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In this paper we study the properties of cold bosons in a two-dimensional optical lattice system 
where Bose-condensation occurs at a momentum point k with non-zero fe-space Berry curvature. 

By combining results from both analytic and numerical approaches, we show that the boson system 
carries non-universal, temperature dependent equilibrium angular momentum and edge current at 
low temperatures. 


It is well known that Berry phase can induce profound 
effect on electronic systeirPl^ and is closely related to 
non- triv ial topological phenomena like the quantum Hall 
effeclP^, quan tum spin Hall effeciP and quantum anoma¬ 
lous Hall effectP^. The topological invariant (Chern num¬ 
ber) characterizing these systems is the integral of Berry 
curvature H(k) over the entire Brillioun zon^. The k- 
space Berry curvature is also predicted to have direct 
consequences on electronic transport ^2121 xhe Berry cur¬ 
vature can be induced by spin-orbit coupling (SOC) and 
several theoretical schemes to realize SOC in cold atom 
systems have been proposed in the last few year^^SHIIl^ 
with some of them being implemented in optical lattice 
experimentalljE^®^. These efforts open a new route for 
scientists to explore plausible Berry curvature-induced 
exotic phases in cold atom systems. For instance, realiza¬ 
tion of t he Ha ldane modeP^, Majorana fermions in ultra¬ 
cold ga^^nHU^ chiral interacting bosonic superfluicP^ 
and the chiral Meissner currents in bosonic ladder^^. 

The combined effect of Bose-condensation and fc-space 
Berry curvature is a natural direction of research in this 
new area. Recently, several related wor ks hav e been done 
within the semiclassical approximatiorP^HSI. in this pa¬ 
per, starting from the Gross-Pitaevskii (GP) equation, 
we study a simple model of weakly interacting bosons 
moving in a lattice with nonzero /c-space Berry curva¬ 
ture at the momentum point where Bose-condensation 
occurs. By combining results from both analytic and nu¬ 
merical approaches, we show that non-universal angular 
momentum and edge current exist in the system at low 
temperatures. The angular momentum has both bulk 
and edge contributions. Furthermore, gapless bulk and 
chiral edge excitations coexist which lead to the modifi¬ 
cation of these entities at non-zero temperatures. 

We start with a continuum model of spinless bosons 
in a two-dimensional optical lattice trapped in an exter¬ 
nal potential P(R). The system is described by the GP 
equatioiP^ = HA! where 

R = -^V? + I4ff(r + Ap(p)), (1) 

yeff(R) = ViR) + UpCR) - fi. 

I4ff(i2) is the effective Hamiltonian acting on the boson 
condensate wavefunction 'k, p(R) = 1I'+(R)1['(R) is the 


boson density and U is the interaction between bosons. 
The optical lattice leads to non-vanishing fc-space Berry 
curvature at the momentum point po (which we shall set 
to be zero in the following) where boson condenses. As 
a result the physical position and momentum (R, p) are 
not equivalent to the conjugate pair (r,p = —iK\/r) any¬ 
more. Instead the physical position is related to r by 
R = r + Ap(p)P where Ap(p) is the fc-space Berry con¬ 
nection and H(p) = Vp x Ap(p) is the Berry curvature 
at the momentum p-point. We note that the existence of 
fc-space Berry curvature at the Bose-condensation point 
implies that Time-Reversal Symmetry is broken in the 
system. 

The kinetic energy term in ([^ should be understood 
as the expansion of the band energy of the bosons in 
optical lattice around band minima po to second order 
in p = —ihVr- To be consistent with this expansion we 
also expand the potential energy term to second order in 
the Berry curvature term Ap, obtaining 

Wff(R) + Ap(p) • Vr + ^(Ap(p) • Vr)"^ I4ff(r). 

The equation can be further simplified if we approx¬ 
imate n(p) H(p = 0), consistent with keeping terms 
up to second order in p in the GP equation. In this case 
we may write Ap(p) = — ^p x $7 (radial gauge), where 
O = n(p = 0)z. We see that Ap(p) is of first order in p 
itself thus to second order in Vr, we only need to keep the 
first order term Ap(p) • Vr in the expansion of I4ff(R)j 
obtaining 

R ^ + ^(A(r) • p + p • A(r)) -k I4ff(r) - p (2a) 

where 

A(r) = -^(f2x VPeff(r)), (2b) 

i.e., the fc-space Berry curvature introduces an effective 
real space vector field A(r) coupled to bosons in the pres¬ 
ence of VI4ff(r) ^ 0. In the following we study the pre¬ 
dictions of this effective GP equation. 

We first distinguish the canonical velocity operators Vc 
and the physical velocity operator v given by 

_ dR _ dr dAp(p) _ dAp(p) 

dt dt dt dt 
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It is straightforward to show that 

Vc = [r, ff] = — + A(r) (3b) 

in m 


and 


V = Vc - ^[p, H] X r2 = Vc + A(r) (3c) 
Zin 

to second order in Vr- 

Correspondingly the real space angular momentum op¬ 
erator is given by 

L = R X (mv) = — -(p X $7)^ x (mv) (4) 

n 

= Lc - —(p • mv) 

where Lc = mr x v is the orbital angular momentum. 
We find that in additional to L^, the presence of non-zero 
Berry curvature leads to an intrinsic angular momentum 
~ f2(p • mv), independent of microscopic details. This 
is consistent with the finding by Chang and NirP, where 
they show that besides center of mass motion, a wave 
packet possesses in general also a self-rotation around its 
center of mass in the presence of A:-space Berry curvature. 
The first and second terms in Eq. Q describes the orbital 
and self-rotation motion of the wave packet, respectively. 

We now apply these results to the boson system. We 
start with the local equilibrium (Thomas-Fermi) approx¬ 
imation (LEA) for the ground state where the solution 
of the GP equation (§ is given approximately by 

^o(r) = VPo(r) = 


for V{v) and 4'o(r) = 0 for /i < P(r). Notice that 
V'eff(r) = V{y) + U p(v) — ^ = 0 at places where 'l'o(r) ^ 0 
and A(r) = 0 has no effect on 4'o(r) in this approxima¬ 
tion. 

To see the effects of /c-space Berry curvature we con¬ 
sider systems with rotational symmetry (i.e. V (r) = 
V(r) where r = |r|) and examining the angular momen¬ 
tum carried by the ground state wave-function. We ob¬ 
tain in LEA 


L = 27ry'rdr4')j'(r)[Lc - y(p • mv)]4'o(r) (6) 


= —Tril 




where p = V(re). Notice that A(r) has no effect on 'l'o(r) 
in LEA and Lc = 0. Asa result there is only self-rotation 
contribution to angular momentum in LEA. 

This expression for angular momentum has an inter¬ 
esting consequence for potential energy of form Vlr) = 
Virjrf) such that Tc measures the size of the system. 
(Notice that the Berry curvature has no effect on 4*0 in 
LEA. As a result, r = R measures the physical positions 


in the system). In this case the condensate wavefunction 
([^ satisfies 

rra 

N = 2Tr rdr\tpo{r)\'^ = {'Kr‘l)n, 

Jo 


where N is the total boson number in the system and n = 
2 Jg xdx\'if!Q{x)\^ is the boson density, x = r/rc- However 
the total angular momentum is given by 


L = 





d'i'oix) 2 

dx ’ 


is independent of the size of the system although N scales 
as the area at fixed density. 

It should be cautioned that the absence of orbital con¬ 
tribution to angular momentum is a special feature of 
LEA which breaks down at the edge of the system. We 
shall see later that in general an orbital edge motion also 
contributes to the angular momentum of the system. The 
edge angular momentum Lc is of order (27rrc) x Tc x 
where jd is the edge current {jd = 0 in LEA). The total 
angular momentum L scales as L ~ a-|-6rc in the presence 
of both self-rotation and edge (orbital) contributions. 

The total angular momentum Eq. suggests that the 
Bose-condensate carries a bulk angular momentum den¬ 
sity 


m(r) = 


n^n d^oir) 2 

2 ^ dr ^ 


( 7 ) 


which can be interpreted as coming from a rotatin g m ass 
current given by j(r) = ^Vr x m(r). Using Eq. M for 
'Lo) we find that LEA predicts that this rotating current 
is given by j(r) = j(r)(/), where (/> is the azimuth angle in 
the two-dimensional plane and 


^ / 1 fdVjr) 2\ 

8mU dr \{p, — V{r)) dr ) 


Notice that j{r) depends strongly on microscopic details 
of the system and has no simple, universal character. The 
rotating current j{r) as a function of i? = r is shown in 
Fig0 with a harmonic trapping potential V{x) = ax'^. 
The current j(r) diverges as i? —t i?c, reflecting the break 
down of LEA at the edge of the boson liquid. Never¬ 
theless, the LEA suggests that circulating current ex¬ 
ists in the ground state of a Boson system which Bose- 
condensate at a momentum point with non-zero fc-space 
Berry curvature. The circulating current concentrates 
mainly at the edge of the sample, and is non-universal 
with structure depending on the details of the micro¬ 
scopic Hamiltonian. 

To go beyond LEA we perform a numerical simulation 
on a tight-binding Hamiltonian with fc-space Berry cur¬ 
vature. We consider interacting spin-1/2 bosons moving 
on a square lattice, 


H = ^ (-t5^_^6R+d,<T + ietz ■ (cr X d)fe^^^6R+d.<T' + h.c.) 
R,,d,fT 

+ Hz ^r_ct^R..o-0’z + UaUfi c- 
R,(T R,,tT 


( 8 ) 
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FIG. 1: The distribution of current from the analytic result 
in a harmonic trap, /r is the chemical potential and is the 
physical boundary of the system. One can see LEA fails at 
the physical boundary. 


where t,a,Bz,U are the hopping energy, Rashba SOC, 
Zeeman energy, the on-site repulsive interaction respec¬ 
tively, and riR = 6^ is the boson density. We apply 

a self-consistent approach by replacing Un^ by U (nR)nR 
where (..) is the ground state expectation value. 

The non-interacting single particle energy spec¬ 
trum is given by (S±(k) = —2t{coskx + cosky) ± 

kx + sin^ ky) + . Expanding around k = 0, 

2 

we see that when B^ > ^, the spectrum has a single 
minimum S_ (k) located at k = 0. In our numerical 
simulation, we consider a weakly interacting boson gas 
in this regime where boson only condenses at the state 
(k = 0) and the ground state is affected solely by the 
Berry curvature at #_(k ~ 0). We consider a finite- 
size, N X N sites sample with open boundary and apply 
the optimal damping algorithm to determine the ground 
state self-consistentljE3. The self-consistent ground state 
boson density distribution and current flow are shown in 
Figj^ We can see that the boson density p{x, y) is nearly 


Density configuration 




FIG. 2: The density configuration (Left) and the current flow 
(Right) in the self-consistent calculation. We take t — 2, 
a = 0.3, Bz = 1, U = 0.2 and = 45 in our calculation. The 
distribution of boson density is consistent with LEA except 
near the edge. The current is chiral and localized around the 
edge. 

uniform in the bulk of the sample and smoothly goes to 
zero at the edge. Indeed the chiral current appears which 
is not found in LEA. This current is mainly localized at 
the edge region where the density p drops rapidly. 

To look at the edge current more carefully we show the 


current jx{y) flowing in x direction as a function of posi¬ 
tion y at fixed x = 23 for several different values of U in 
Figj^ We see that the edge region becomes smaller and 
sharper when the interacting parameter U increases. We 
emphasize here that the chiral edge current we observed 
is qualitatively different from the ones in gapped topo¬ 
logical systems :(I) the Chern number in our system is 
not required to be nonzero, since the edge current is re¬ 
lated to Berry curvature at a local /c-point only; (2) our 
system is a gappless system, owing to the spontaneous 
U{1) symmetry breaking. 



FIG. 3: The charge current distribution in our model Eq.Q 
for several values of interaction strength U. The value of other 
parameters are the same as in Figj^ The current is mainly 
localized at the edge and the localization is strengthened by 
increasing U. 

To examine the scaling relation between angular mo¬ 
mentum and sample size A ^ i?c, we fix the average 
density n = 2 per site in our simulation. The L — 
relation is shown in Fig|^ for several values of U. We 
see from Fig|^ that the total angular momentum scales 
L ~ a + bR^, consistent with our analytical result. Notice 
b depends on 17, consistent with our expectation that L 
is non-universal and depends on the microscopic Hamil¬ 
tonian. 



FIG. 4: The scaling relation between L and size dimension 
A ~ i?c for different interacting strength. We find that the 
total angular momentum scales as L ~ a -I- bR^. Notice 6 is a 
function of U. 

Next we consider the effect of Berry curvature on the 
low energy excitations. For trapped bosons, there are 
bulk excitation modes and excitation modes localized 
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near the edge of the boson cloud . The excitatio ns can 

be obtained by writing 4'(r) = ^/po(r) + 5p(v, 

and solves the linearized GP equatiorP^ for 6p and 69. 

Keeping terms to order V^, we obtain after some 
straightforward algebra the well known hydrodynamic 
equation^ 

^ + V • (poVc) = 0 (9a) 

ot 

and 

m^=-^{V{v) + U5p)-C^W{^)xn) (9b) 

where Vc = —V56 + A = dr/dt, in agreement with Eq. 

At the bulk of the sample, W = 0, po is a constant 
and it is straightforward to obtain the usual Goldstone 
mode Sp{r,t) = with dispersion 



m 


Notice that the Berry curvature has no effect on the bulk 
dispersion. 

The situation is rather different in the surface region. 
By approximating the potential in the surface region (in 
x-direction) as a linear function of the coordinate V (r) — 
/i ~ Ex and assuming the LEA for po, we obtain for 
q « where 6 ~ {h^/2mF)^, 

Sp{x, y, t) = /( 9 x)e-«l"l±*«^e-“‘.(<z > 0) (10) 

and ±g represents two different modes propagating along 
positive and negative y direction, respectively. 

The solution for Eq. (|^ is the Laguerre polynomials 
Ln{2qx) with dispersion relation 


'■W) = 


771 16 4 


( 11 ) 


consistent with the result obtained by Price and 
Coopeil^. Note that for id 7 ^ 0, the excitation energies 
correspond to the left and right propagating modes are 
different. 

The surface excitation results in a T 7 ^ 0 contribution 
to the equilibrium surface current given by 


jiT) = (nB(w+(g))z;+(g)-nB(w„(g))u„(g)) 

n,q>0 

( 12 ) 

where v^{q) = WqUj^{q) is the group velocity and 
is the Bose-Einstein distribution function. Using the fact 

that ?;n^~^((z) = in ID; we obtain 

27rr / /"“a /"“a 

js{T) = f dw+( 9 ) 713 ( 0 ;+) - duj-{q)nB{ujq) 

~ rc(fc3r) , (13) 


where js{T) depends on the high momentum cutoff. Here 
we only consider n = 0 contribution but similar conclu¬ 
sion can be drawn for n 7 ^ 0 . 

The direction of this finite temperature current is op¬ 
posite to the zero temperature edge current, meaning 
that the equilibrium edge current is reduced as temper¬ 
ature increases. On the other hand, the total angular 
momentum carried by the system is a sum of two parts, 
the bulk (self-rotation) contribution Lb{T) and the sur¬ 
face (orbital) contribution ALs{T) ~ Tc x js{T). The 
bulk contribution at low temperature is given by 

o -jrA 

ALb{T) = -2 ^ <fq{hqfnB{uJq) (14) 

^rlm?{kBT)'^ 

^ mplu^ ’ 


where a = 


0 e“-l 


dx. 


In particular, we find from Eqs. (131 and ( |14[ ) that the 
low temperature correction to total angular momentum 
AL(T) = L(T)-L(0) = ALb(T)+ALs(T} is dominated 
by the bulk (self-rotation) contribution at low tempera¬ 
ture and total angular momentum of the system increases 
as function of temperature at low temperature although 
the edge current decreases as temperature raises. 

Experimental Implementation To realize our model 
experimentally, we need to ( 1 ) create pseudo-spin- 1/2 
bosons, (2) implement Rashba SOC, and (3) generate 
Zeeman splitting in the pseudo-spin space. The pseudo¬ 
spin 1/2 bosons can be created as in the experimental 
setup of Spielman’s grouffi^ which breaks the degener¬ 
acy of ®^Rb, E = 1 ground state manifold and select 
the two internal spin states with lowest energy to be the 
basis of pseudo-spin-1/2 particles. Secondly, to generate 
2D Rashba SOC, one may employ the proposal of Liu et 
to apply two Raman fields induced by four lasers. 


The two Raman fields generate (sin/coa; -t- i sin fco 7 /)o|a^ 
term in the Hamiltonian which is 2D SOC. Thirdly, one 
needs to generate a large Zeeman energy splitting in this 
pseudospin basis. The small two-photon off resonance 
term mentioned in Ref.^ may not be sufficient and an¬ 
other counter propagating Raman field is needed to real¬ 
ize Rabi oscillation as in Spielman’s experimen!^. Eesh- 
bach resonance can be applied to control the boson-boson 
interaction by tuning s-wave scattering length as a func¬ 
tion of external magnetic fielcP^. 

Last but not least, to detect the edge current in our 
system, one may conduct the time-of-flight(TOF) mea¬ 
surement for two spin species separated by Stern-Gerlach 
effeclP^. With SOC and interaction, we find a cross-like 
distribution around T point when analyzing the density 
of minor component 77 ^(k) (shown in Eigj^. The maxima 
at finite k points are results of the edge mode. The chi- 
jral edge excitations can be detected directly by the Bragg 
'scattering process which based on a two-photon process 
that directly transfers energy and momentum to an en¬ 
semble of atomj ^. Ano ther possibility is to use quantum 
quench techniqued^^^ which detects the structure of the 
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FIG. 5: The density distribution of the minor component in 
momentum space n^(k). Without SOC the density of the 
minor component at ground state is almost zero due to the 
Zeeman field splitting. With SOC and interaction U = 0.2, 
the density distribution of n|(k) spreads out to higher mo¬ 
mentum, with four maxima at finite k points corresponding 
to edge mode. 

boson wave function. One may apply a shaping poten¬ 
tial well which is initially prepared to have a surrounding 
edge current. After removing this shaping potential, one 


can measure the cloud density revolution, revealing the 
information of angular momentuirP^. 

Summarizing, we study in this paper the Gross- 
Pitaevskii (GP) equation for a simple model of weakly 
interacting bosons moving in a lattice with nonzero k- 
space Berry curvature at the momentum point where 
Bose-condensation occurs. By combining results from 
Local Equilibrium (Thomas Fermi) approximation and 
numerical self-consistent approach, we show that non- 
universal angular momentum and edge current exist in 
the system at zero temperature. The angular momentum 
has both orbital and self-rotation components. Further¬ 
more, the equilibrium edge current is reduced by gapless 
chiral edge excitation as temperature increases but the 
total angular momentum carried by the system increases 
as temperature increases as a result of bulk self-rotation. 
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by HKRGG through grant No. HKUST/CEF/13G. 
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